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Microresonator based optical frequency combs (MFC) have demonstrated promise in extending the 
capabilities of optical frequency combs. Here we demonstrate all optical stabilization of a low noise 
temporal soliton based MFC in a crystalline resonator via a new technique to control the repetition 
rate. This is accomplished by thermally heating the microresonator with an additional probe laser 
coupled to an auxiliary optical resonator mode. The offset frequency is controlled by stabilization 
of the pump laser frequency to a reference optical frequency comb. We analyze the stabilization 
by performing an out of loop comparison and measure the overlapping Allan deviation. This all 
optical stabilization technique can prove useful as a low added noise actuator for self-referenced 
microresonator frequency combs. 


Since their inception |T], optical frequency combs 
(OFC) have significantly changed the way precise optical 
frequency measurements are performed. Typically cre¬ 
ated from a train of ultra short pulses emitted by mode- 
locked lasers in combination with supercontinuum gener¬ 
ation an equidistant grid of optical frequency markers is 
generated. Any of the frequency components can be writ¬ 
ten in terms of the expression f n = n- / rep + / ceo . Where 
n is an integer and parameters / rep and / ceo correspond 
to the repetition rate of the pulse train and to the carrier 
envelope offset frequency. Discovered in 2007, microres¬ 
onator based optical frequency combs (MFCs) [2 [3] gen¬ 
erated via parametric frequency conversion HUS] opened 
up a new parameter space for OFCs, achieving compact 
form factor and on chip integration, enabling high repeti¬ 
tion rates (typically > 10 GHz), and broad spectral from 
the near to the mid-infrared. Such large mode spacing in 
the microwave regime [3] is beneficial in applications such 
as astronomical spectrometer calibration |6j, dual comb 
coherent Raman imaging [7], high speed optical sampling 
and optical telecommunication [5 . 

A detailed understanding of how frequency comb for¬ 
mation arises in these systems has been obtained in paral¬ 
lel by both experimental [9j and theoretical work puma. 
To generate a MFC, a continuous wave laser is tuned into 
a cavity resonance and once the intracavity intensity is 
above the parametric threshold frequency comb forma¬ 
tion begins via degenerate and non degenerate four-wave 
mixing and many possible states can result. Multiple 
sub-combs are formed when the variation in the cav¬ 
ity free spectral range is small compared to the cavity 
linewidth. These states exhibit many types of noise [9], 
but can be brought into the low noise regime via S — A 
matching |5J, parametric seeding m, injection locking 
[IS HU, or via the observation of mode-locking [15]. In 
the opposite regime it is possible to form low phase noise 
states for narrow bandwidth combs mmm- An addi¬ 
tional route that is used in this work is temporal dissipa¬ 
tive cavity soliton formation PHD!. Temporal dissipa- 
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tive soliton formation results in deterministic low noise, 
smooth spectral envelope frequency combs that in time 
domain give rise to ultrashort optical pulses. In this man¬ 
ner sub 200 fs pulses were generated in a crystalline mi¬ 
croresonator using a single soliton state, and in conjunc¬ 
tion with external broadening it was possible to deter¬ 
mine / C eo via a 2f-3f self-referencing m- 

The ability to self-reference an OFC is one of the key 
requirements of a MFC being useful for precision metrol¬ 
ogy. In some applications of MFCs such as atomic clocks 
m it is important to also stabilize both / rep and / ceo . 
It has been well established m that it is possible to 
independently control both parameters. One unique as¬ 
pect of MFCs is that the pump laser (/ Q ) is the central 
comb line and changing its frequency changes / ceo . This 
can also affect the repetition rate via thermal and Kerr 
effects but an orthogonal basis can be found with other 
control parameters [22]. Controlling the repetition rate is 
possible through a variety of means: changing the pump 
power m, actuating on a piezo electric crystal that is 
in contact with the resonator [ 23 :, or heating or cooling 
the whole system [23]. Here we demonstrate all optical 
stabilization of both / rep and / ceo of a soliton based MFC 
by controlling f Q via the pump laser detuning and a new 
optical technique for the control of / rep . For this, an 
additional probe laser is coupled into a different cavity 
mode than the ones used for frequency comb generation. 
By adjusting the detuning of the probe laser, the power 
coupled into the microresonator is changed affecting / rep 
(see figure [l]) . Auxiliary modes have also been used to 
monitor the resonator temperature and thus the repeti¬ 
tion rate (25] and to compensate for thermal nonlinear 
effects [26]. 

Whispering gallery mode microresonators (WGM) 
support azimuthally symmetric optical modes, and in 
this work the modes are confined in a small protrusion 
around an axially symmetric crystalline magnesium flu¬ 
oride (MgF 2 ) microresonator [T5] 27] , The geometry of 
the resonator, along with the material properties deter¬ 
mine the number and type of optical modes that can 
be supported. It is possible to engineer the protrusion 
to support only one mode [28]. However in this work 
we leverage having more than one spatial mode in the 
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Figure 1. Microresonator and actuation scheme: a. The 
MgF 2 WGM microresonator is formed by the top protrusion 
going around the circumference (the lower one is not used 
in this work). A tapered optical fiber couples in and out the 
pump and probe lasers and a portion of the generated soliton. 
b. Schematic of the basic concept used to stabilize the repeti¬ 
tion rate. The CW pump laser (blue) is tuned into resonance 
of the high Q modes that produces the soliton frequency comb 
producing mode (blue). Tuning the pump frequency changes 
f 0 . The auxiliary lower-Q spatial mode family of the cavity 
is shown in red. The probe laser is tuned to the high fre¬ 
quency side of one of the resonances, c. A zoom in of the 
relevant auxiliary mode. Tuning the probe laser in the res¬ 
onance changes the power coupled in and thus the FSR for 
both modes and changes / rep . 


resonator for stabilization of the repetition rate. One 
advantage of crystalline resonators is the exceptionally 
high quality factors (Q) exceeding 10 10 [29) [30] that can 
be obtained. In this work the resonance used for OFC 
generation has a Q ~ 10 9 and a free spectral range of 
~ 14.09 GHz. To form the OFC, continuous wave laser 
light from a 1553 nm fiber laser (pump laser) with 240 
mW is coupled evanescently into the resonator via a ta¬ 
pered optical fiber m- Single or multiple soliton states 
are generated in the resonator by performing a controlled 
frequency scan in the direction of decreasing frequency 
with respect to the resonance , and by choosing the cor¬ 
rect stopping frequency different soliton states can be 
accessed [IB . A portion of the light is coupled out via 
the tapered optical fiber, and the resulting spectrum for 
a single soliton can be seen in figure [2j 

The principle of the stabilization scheme for / rep and f Q 
is presented in figure [l] f Q is stabilized by offset locking 
the pump laser to an external reference [22 . Actuating 
on / rep is done via a second “probe” laser at ~ 1552.1 nm 
with ~ 2 mW, coupled and locked to a resonance of an¬ 
other mode family with lower Q. By changing the lock 
point in this resonance, more or less optical power is cou¬ 
pled into the resonance. Due to absorption losses in the 
cavity this heats the resonator, changing the overall size 
and refractive index and thus the repetition rate with¬ 
out loss of the soliton state. In addition the Kerr effect 
also changes the effective index of refraction on fast time 


scales. However, the magnitude of this effect is negligible 
compared to the thermo-refractive index change. 

The relevant stages of the experimental setup for the 
stabilization of the soliton frequency comb are presented 
in figure [2j The locking setup for f Q is shown in the 
green panel. A portion of the pump light is heterodyned 
with light from a reference comb consisting of a com¬ 
mercial self-referenced and stabilized fiber based optical 
frequency comb with a repetition rate of ~ 250 MHz. 
An optical bandpass filter (OBP) is used to filter (~ 50 
GHz bandwidth) the light preventing saturation of the 
photodetector (PD). The heterodyning produces multi¬ 
ple beat signals between the pump laser and the fiber 
frequency comb, which represent the frequency offset of 
the pump with the near by comb lines. One of the beat 
signals is electronically mixed to 20 MHz with a local 
oscillator phased locked to the master oscillator (com¬ 
mercial atomic clock). A portion of the signal is sent to 
a frequency counter. The rest is sent to a digital phase 
comparator (PC) with a 20 MHz reference, which outputs 
an error signal to a proportional and integral (PI) servo 
controller (servo 1). This provides feedback to change 
frequency of the pump fiber laser via an internal piezo¬ 
electric actuator. 

The stabilization setup for / rep is shown in the orange 
box in figure [2^i. The soliton pulse leaving the resonator 
sits on a strong continuous wave background, which is 
partially filtered out using an optical notch filter (ONF 
- 50 GHz bandwidth) to prevent saturation of the sub¬ 
sequent high speed PD. To control / rep , the additional 
probe fiber laser is coupled into the same tapered opti¬ 
cal fiber as the pump laser, and its frequency is tuned 
into a different spatial optical mode resonance of the mi¬ 
croresonator. A portion of the light coming out of the 
resonator is sent to the probe locking setup depicted in 
the blue box in figure [2^i. An OBP filter is used to pass 
primarily the transmitted probe light to a photodetec¬ 
tor. The probe laser’s frequency is then locked to the 
high frequency side of the cavity resonance, where ther¬ 
mal locking is also supported. To suppress nonlinearities, 
lower powers and a lower Q mode are used. It was not 
necessary to determine the exact mode of the cavity, be¬ 
cause many different modes showed control over / rep . A 
PI servo controller with an adjustable set point (servo 2) 
controls the lock point of the probe laser on the side of 
the cavity resonance. The servo controller feeds back to 
a piezoelectric actuator on the probe fiber laser, which 
adjusts the laser’s frequency. To fix / rep an additional 
control signal is needed which is derived from the repeti¬ 
tion rate heterodyne signal described above. This signal 
is down mixed from 14.094 GHz to 20 MHz using a local 
oscillator that is referenced to the master oscillator. The 
signal is filtered with an electronic bandpass filter and a 
portion is sent to a counter. The rest is sent to a digital 
PC with the 20 MHz signal from the master oscillator 
where any phase error produces an error signal for a PI 
servo controller (servo 3). The correction signal from the 
output of servo 3 is input to the servo 2, which maintains 
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Figure 2. All optical stabilization of a soliton frequency comb. a. Experimental setup for generating and stabilizing the 
soliton frequency comb. The details are provided in the text. The components are: continuous wave (CW) pump and probe 
lasers, erbium doped fiber amplifier (EDFA), fiber polarization controller (FPC), optical bandpass filter (OBP), optical notch 
filter (ONF), optical spectrum analyzer (OSA), photodetector (PD), local oscillator (LO), phase comparator (PC), electrical 
spectrum analyzer (ESA) b. The generated optical spectrum as measured on the OSA. c. A zoom in of the optical spectrum. 
The pump laser and the probe laser’s position is shown. Due to the limited resolution of the OSA the probe appears merged 
with the 7th comb line. The position of the out-of-loop optical filter, covering several comb lines of both frequency combs, is 
shown overlaid on the spectrum. It is centered — 73 comb lines away from the pump laser (Oth comb line). 


the probe laser’s lock on the probe resonance, adjusting 
the lock point detuning of the probe laser to maintain 
the repetition rate phase lock. The time constants of the 
systems were not directly measured but rather the ap¬ 
propriate proportional and integral time constants were 
experimentally determined. It should be noted that a sub 
1 s integral time constant was used for the / rep servo con¬ 
troller due to the observed slow response of the system, 
indicating that thermal effects appear to dominate the 
response. To verify the stabilization is not injecting sig¬ 
nificant noise an independent out-of-loop measurement 
is performed by heterodyning a portion of the generated 
soliton comb with the reference fiber based comb (gray 
box in figure [2J. An OBP filter is used to select a por¬ 
tion of both frequency combs centered at 1544.6 nm. Due 
to the — 50 GHz bandwidth width of the filter multiple 
heterodyne beat signals are observed. Using a tunable 
electronic filter, one signal f Q \ is selected and sent to a 
frequency counter. 

We analyze the stability of the full system by perform¬ 
ing Allan deviation measurements on / rep , / Q , and f Q \ on 
both a single and multiple soliton states. The hetero¬ 
dyne beat frequencies are measured on Hewlett Packard 
53131A high resolution counters and the overlapping Al¬ 
lan deviation (OAD) is processed from the recorded fre¬ 
quency series. For the single soliton state shown in fig¬ 
ure IT counting data was taken for 3086 s at 100 ms gate 
time, yielding the OAD plot in figure [3j At 100 ms the 
absolute fluctuations of the 14.094 GHz rep rate beat are 
<TA,rep = 4.98 Hz corresponding to a fractional deviation 
of 3.5 x 10 10 . The absolute frequency fluctuations of 
f Q at 100 ms are cta, 0 = 8.9 x 10“ 1 Hz. Only the abso¬ 
lute fluctuations of f Q can be measured since the exact 
value of / 0 is not known in this experiment. For out- 
of-loop measurement the beat frequency was 70.3 MHz 


and the absolute fluctuations at 100 ms were found to be 
<ta,o 1 = 283 Hz. The bump in <JA 5 rep and cr a, 0 i is probably 
a result of the slow thermal response of the system. Per¬ 
forming a fit to the slopes, using the data for times > Is 
to make sure the system was stabilized, we find that / rep , 
/o and /oi average down like T~p- 53 , t” 0,48 and r^ 0,53 . 
We observe here a dependence of ^ r -0,5 while for a 
phase locked system the expected dependence is rather 
~ t - 1 [32]. This observation is related to the properties 
of the Hewlett Packard 53131A counters used. They are a 
A-type counter because they perform a weighted average 
of the frequency over the gate time to enhance the resolu¬ 
tion. As a result, the computed deviation differs from the 
true OAD but still measures the stability of the system 
m ■ Also, this counter cannot be read without dead time 
between consecutive measurements. This creates a bias 
]341 in the processed OAD for the phase-locked system 
and leads to a r -0,5 dependence in the presence of white 
phase noise. In our case the decreasing behavior of the 
fluctuations for longer gate time still shows the stability 
transfer of the master oscillator to the soliton frequency 
comb. The out-of-loop measurement is approximately a 
factor ~ 71.5 higher than a\,rep for times > Is. The out- 
of-loop OBP filter is centered approximately n Q \ ~ 73 
comb lines away from the central pump line. We can 
write = An 0 i x / r ^ p + /£ with k being the measure¬ 
ment number and where A n Q i is the mode number rel¬ 
ative to the central comb line, and thus a\ Q i = G \ 0 + 

An ol°A,rep + An °l {(/rep 1 " frep)(fo + 1 ~ /o j) and since 

cr A,o < cr a, rep for r > 1 s we expect cr A , 0 l ~ An 0 i cr A , r ep 
which in good agreement with the observed offset. 

We also demonstrate it is possible to fully stabilize 
multi-soliton states. The desired state is created using 
the same tuning technique m for the single soliton. 
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Figure 3. Overlapping Allan deviation for a stabilized single 
soliton state. Counting data with 100 ms gate time was taken 
for 3086 s. The resulting OAD is shown by the traces for the 
pump frequency (blue), repetition rate (red), and the inde¬ 
pendent out-of-loop measurement (yellow). The error bars 
are one standard deviation. The dashed lines are a linear fit 
for the data >1 s. The gray trace shows the result of counting 
the 20 MHz master oscillator used in the experiment. 


The spectrum for the created state is shown in figure [4] 
Having multiple solitons alters the microresonator optical 
spectrum as a result of interference between the different 
frequency components. The time domain pattern of the 
intracavity amplitude can be reconstructed by fitting an 
analytical expression to the optical spectrum [55] , The 
left inset in figure [4] shows only the peak amplitudes for 
clarity. It is estimated that the state in figure [4] is a 
two soliton state spaced apart by ~ 90° in the microres¬ 
onator. The same stabilization was applied as for the 
single soliton with ~ 3000 s of counting data taken at 1 s 
gate times and the OAD was calculated for / rep , f Q , and 
/oi, shown in figure [4|d. At 1 s, the absolute fluctuations 
of rep rate beat is cr a, rep — 14 Hz corresponding to a frac¬ 
tional deviation of 7.8 x 1C) -11 . The absolute frequency 
fluctuations of f Q at 1 s are <ja, 0 = 5 x HR 2 Hz. For 
the /oi the beat frequency was 48.7 MHz and the abso¬ 
lute fluctuations at 1 s were found to be cta, 0 i = 79.5 Hz. 
This is a factor of ~ 72.6 higher than fluctuation seen in 
/ rep , which again agrees well with the out-of-loop filter 
position. The bump in <7A, r ep and cta, 0 i is absent here 
as the gate time was taken to be 1 s. Fitting the data 
gives the slopes r~ p ,5 °, t“ 0,48 and r^ 0 ' 50 . This data has 
a lower initial OAD. The servo parameters were changed 
between the single and multi-soliton case, which poten¬ 
tially explains the difference. 

In summary we have demonstrated that an auxiliary 
laser can be used as actuator to stabilize the optical fre¬ 
quency comb spectra associated with a temporal optical 


soliton in a microresonator. The distinct advantage of 
this actuator is that it is easily implemented, and that 
in addition the added noise is only limited to laser noise. 
Our technique is applicable to a wide range of other mi- 
croresonator platforms. 
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Figure 4. Multi-soliton spectrum and overlapping Allan de¬ 
viation. a. A measurement of the optical spectrum from the 
OS A of the stabilized two soliton state (blue). The peaks 
of the spectrum are detected and fitted with a theoretical ex¬ 
pression [35] (left inset). This yields the duration, the number 
and relative position of the solitons in the cavity (see the inset 
in the upper right), b. Counting data at a 1 s gate time was 
taken for 3000 s. The resulting OAD is shown by the traces for 
the pump frequency (blue), repetition rate (red), and the in¬ 
dependent out-of-loop measurement (yellow). The error bars 
are one standard deviation. The dashed lines are a linear fit 
to the data. The gray trace shows the result of counting the 
20 MHz master oscillator used in the experiment. 
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